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Acronyms 


1D One Dimensional 

2D Two Dimensional 

3D Three Dimensional 

a Solar Absorptance 

Alt Altitude 

ASSY Assembly 

Bij Interchange Factor 

BOL Beginning of Life 

CAD Computer Aided Design 
Cp Specific Heat 

DT Temperature Difference 
E Infrared Emissivity 

6* Effective Emissivity 

EOL End of Life 

FD Finite Difference 

FE Finite Elements 


GOLD Goddard Open Learning Design 

GSFC Goddard Space Flight Center 

Ncona/Meyap Condensor / Evaporator Heat Transfer Coefficient 
HP Heat Pipe 

HR Heat Rate 

ID Identifier 


IF 
Incl 
IR 


Interface 

Inclination 

Infra Red 

Thermal Conductivity 

Monte Carlo Ray Trace 

Multi-Layer Insulation 

Non-Condensable Gas 

Phase Change Material 
Proportional-Integral-Derivative 

Density 

Right Ascension of Ascending Node 
Radaition coupling 

Reservoir 

Spacecraft 

Systems Improved Numerical Differences Analyzer 
Standard for the Exchange of Product Model Data 
Structural-Thermal-Optical-Performance 
Thermal Desktop 

Thermo Electric Cooler 
Ultraviolet-Visible 

Variable Conductance Heat Pipe 

View Port 


PURPOSE 


This presentation is meant to be an overview of the model building process 


It is based on typical techniques (Monte Carlo Ray Tracing for radiation exchange, Lumped Parameter, Finite 
Difference for thermal solution) used by the aerospace industry 


This is not intended to be a “How to Use Thermal Desktop” course. It is intended to be a “How to Build Thermal 
Models” course and the techniques will be demonstrated using the capabilities of ThermalDesktop. Other codes 
may or may not have similar capabilities... 
The General Model Building Process can be broken into four top level steps: 

1. Build Model 

2. Check Model 


3. Execute Model 


4. Verify Results 


BUILD MODEL: General Model Building Process 


1. Define Thermo-Optical and Thermophysical Properties 
2. Get CAD from Mechanical Designer — STEP format preferred, make sure it is solids and not only surfaces 
Process CAD into layers for controlling visibility 
3. Build and define each CAD component 
Geometric Surfaces, Nodes, Elements (Plate, Solid) - Add Comments to help 
For Solid Element geometry, apply surface coat of zero thickness plates 
Determine Node Subdivisions and Edge vs. Centroid Nodalization 
Merge Nodes 
Assign Node Numbers and Submodels 
Assign Material and Thickness (Orienters for Solid Elements) 
Assign Active Sides, Optical Properties, and Radiation Analysis Groups 
Assign MLI 
4. Determine Variable Parameters 
Add Symbols (Thickness, Power, Multipliers, IF Values, etc) and assign to geometry if needed 
5. Accommodate Moving Geometry 
Define Assemblies and/or Trackers and Assign Geometry (Include nodes to move FE) 
6. Establish Model Connectivity 
Add Contactors between surfaces/edges at interfaces 
Add Conductors between nodes/surfaces 
7. Establish Boundary Conditions 
Add and Define Heat Loads 
Add and Define Heaters 


BUILD MODEL: Define ThermoPhysical and ThermoOptical Properties 


Edit Optical Properties 


Current Optical Property Database: 
version 1\thermal model from 2015 for hume\thermal desktop model\rcoptics.rco 


¢ Add ThermoOptical Properties: 


New property to add: | Add 
Name he Solar Absorptivity — IR Emissivity a/e 
O.1emiss 1.000 0.100 10.000 
AgFep 0.130 0.740 0.176 
Black_Anodize 1.000 0.780 1.282 
BlackPaintZ306 0.940 0.900 1.044 
GBK 0.640 0.740 0.865 
M55) 0.930 0.740 1.257 
Solar_Cell 0.910 0.810 1.123 
293_White_Paint 0.160 0.920 0.174 
2Z93P_WhitePaint 0.200 0.880 0.227 
| Import 
Gel Geo 
Edit Optical Property - BlackPaintZ306 
Comment: = | . Set Color... | 
Use Properties: | Basic Props for Radks and Heat Rate Calculations NA } 
Basic | \Wavelenath Dependent 
Solar 
Absorptivity: 0.94 Vs. Angle Vs. Temperature 
Transmissivity: 0 Edit Table... Vs. Angle 
Specularity: 0 Edit Table... } Vs. Angle 
Transmissive Specularity: 0 Edit Table... | Vs. Angle 
Refractive Indices Ratio: 1 
Infrared 
Emissivity: 09 Edit T able Vs. Angle Vs. Temperature 
Transmissivity: 0 " Edit Table =a] Vs. Angle 
Specularity: 0 Edit Table... } Vs. Angle 
Transmissive Specularity: 0 Edit Table... | Vs. Angle 
Refractive Indices Ratio: 1 
| Cancel Help 
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Open/Create Property Database 

Add Property and Specify Name 

Define a andé 

Include comment/source 

Repeat until all properties added... 

Repeat for EOL or other cases 

Can import from other models 

Define any Property Aliases that may be used 
(e.g. Radiator_Coating). Use ALIAS_ prefix 


Add ThermoPhysical Properties: 


v 
v 
v 


Open/Create Property Database 
Add Property and Specify Name 
Define p, C,, and k (or k,, k,, k,) 

Y For MLI, define s* 

Y For PCM, define FUSION 
Can import from other models 
For MLI, recommend using p of 0.6 kg/m? and 
k and C, of 0. For MLI assignments, define 
thickness for MLI as 1.0 m 

VY This will allow mass of MLI to be 

estimated based on 0.6 kg/m? 

Define any Property Aliases that may be used 
(e.g. Panel_ Material). Use ALIAS__ prefix 


wn 


Edit Thermophysical Properties 


Current Thermophysical Property Database: 


\wersion 1\thermal model from 2015 for hume'thermal desktop model\tdthermo.tdp 


New property to add: | ‘Add | 
Name . Cond [W/m/C] Dens [kg/m*3] Cp L/kg/C] Eff Emiss Type 
416061-T6 167 2700 896 
Al 7075-17351 155 2810 360 
AL_6061_T6_Core_3.1pcf 1.9685 48.5041 300 
APG 668 2700 896 
Copper 385,827 8960 385 
e*nom it) 06 0 0.02 
G10 0.29 1300 600 
K1100 500 1 1 
M55) 36 1633.11 800 
Stainless Steel 16 8030 504 
VT4?_board 05 1900 600 
{ Edit | Delete | [ Copy [ Rename Import 
Cancel [Help 
Edit Thermophysical Property - Al 6061-T6 
Property: Al 6061-76 | Set Color... 
Comment: * 
Use Properties: {Basic Properties for Material x] 
Basic 
Conductivity [W/m/C] 
k 167 Edit Table... Use Table Pressure... Use Pressure Scale: 1 
@ Isotropic 
Anisotropic 
SpecificHeat [J/ka/C] 
@ 896 Use Table Fusion... Use Fusion 
Density [kg/m*3] 
tho 2700 Scale: 1 
Effective emissivi ity 
estar 0 (used for insulation and core) 
Recession 
Allow Recession 
Cancel} Help 
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BUILD MODEL: Importing CAD 


(A Import File 


¢ To import a STEP file: 

Y Type import, select file, and wait for notification in bottom right 
(Known issue with this command with AutoCAD 2014 in 
Windows 10) 

VY Type explode and select object (this will explode the block into 

constituent parts) 

VY Type explode followed by ‘filter (‘makes this command transparent 

within explode command) FENGINEER Eee 

V Select Block in the Select Filter DropDown and Click Add to Listif guzeueme Saat 

not already there ™ _ ee 

Y Click Apply and type all at command prompt 

Y If number of objects found is not zero, hit return twice (to end filter I — 

and to end explode) a 

Y Repeat the Explode/Filter combination until zero objects are found | 

Y At this point, everything should be a 3DSOLID or non-block entity 

type 
Y Create layers for all the components for the design (/ayer) fo Seren ayer SexchiorlyerQ 


Look in: [ ju) Steinfela >| © &@ MH vews ~ Too wv 


) = . Name " Date modified Type 
|_| Sensor_block_v2.STEP 8/29/2017 12:03 PM _— STEP File 


3D Studio (°.3ds) 
AC 


CATIA V. 4 (" model:*.session;*.exp;*.div3) 
CATIA V5 (°.CATPart;*.CATProduct) 


Select Filter 


Edit tem | | Delete Clear List 


. Bl & & | <3 EK GHP 

Y All Solids can now be placed on layers and have colors changed to 2 ites BO ee ct gre Sl hee eee ed 
. Be ee ay all res Ww 0 y ‘ af ws Continu... — Defa... 0 Colac =) ms 
help with visibility control nN rr ng hr TS 


Y Solids can also be copied and sliced (slice) to take measurements 
and “see” what is inside 
VY Solids can have a material applied to them make them transparent 


Invert filter «|| « 


(th) Layer Properties Manager 


All: 3 layers displayed of 3 total layers 


BUILD MODEL: Working with CAD ‘ 


¢ Many functions help with using CAD: 


MA Se A SSS 


ZOOM: enlarge and shrink image on screen 
3DFORBIT: rotate image in screen 

OSNAP: allow user to select midpoint, endpoint, quadrant, intersection, node, center, perpendicular, nearest 
UCS: define 2D working plane 

MOVE: translate geometry 

COPY: copy geometry within drawing 

ROTATE3D: rotate geometry about specified axis by specified angle 

MIRROR3D: create object that is mirror image of selected object about specified plane 
3DALIGN: align geometry based on 3 points in reference frame and destination frame 
SLICE: cut solid objects by slicing plane 

LINE: create line segment 

DIVIDE: create N points along line segment 

EXPLODE: convert assembly into selectable constituent parts 

DIST: measure distance and angle between two points 

ARRAYCLASSIC: create copies or an object as array (rectangular or polar) 
PROPERTIES: change layer, color and linetype for selected objects 
LAYER: allow visibility, color, linetype changes for specific layers 

ORTHO: toggle if second point is orthogonal (along X or Y) to first point 
FILTER: allow removal or inclusion of selection objects based on properties 
FILLET: connect two, non parallel lines 

ERASE: remove object from drawing database 


BUILD MODEL: Building eaanenns 


Most models represent 3D CAD with 2D simplifications 

The representations include additional properties that make them thermal | 

— Thermo-optical properties and radiatively active sides 

— Thermo-physical properties and thickness 

— Node assignments and subdivision : 

— The node assignments form the basis for a network thermal model which 
is used to compute temperatures and heater powers. It is the link between 
the geometric and thermal models 

Some surface shapes follow naturally occurring geometry 

— Rectangle, Cylinder, Disk, Cone, Triangle, Sohere, Paraboloid, etc. 

Other 2D shapes do not conform to simple geometry 

— For these, the use of finite elements is recommended, which 
allows a user to more closely match CAD geometry 

Sometimes, it is necessary to model heat flows as three-dimensional 

— Solid Primitives (Box, Disc/Cylinder, Sohere, Cone) exist which provide the 
same properties as their 2D counterparts, but allow for a network based on 
3D subdivision 

— For more complex shapes, 3D finite elements (Tetrahedron, Wedge, Brick, 
and Pyramid) are available. Note that these do not include radiation 
properties... 


BUILD MODEL: Meshing Geometry =e 


e Use of Edge Nodes recommended over centroid nodes 

Y Centroid and edge nodes both account for conduction within a single surface 

vY Edge nodes move nodes to edges of each sub-surface instead of centroid 

vY Edge nodes that are merged allow conduction between surfaces via the 
shared edge nodes 

VY Use of finite elements allows for more complex shapes to be modeled 

Y Building elements by hand not as time consuming as may be thought 
> Define nodes first by placing user nodes at location on CAD geometry 
> Define element by selecting nodes in right hand rule order to define “top 

side” 

Y Solid elements can also be defined and allow non-isotropic properties 
> Manually define by selecting nodes (one face then the opposite) 
> OR Revolve/Extrude Plate elements 

Y Benefit of elements is ability to move nodes to desired locations 

Y Elements cannot be moved... 
Elements defined by nodes... 
Move the nodes instead 

VY Shell coat solid elements 

Y Define surfaces as either edge 
Node or centroid node 

v Merge nodes 

Y Assign rest of properties... 

* Extremely important!!!! 


BUILD MODEL: Working with Finite Elements 


Finite Elements are a powerful tool, but can be more difficult to work with than primitive surfaces 

Elements cannot be double sided and the node numbers are not able to be changed by the element. You can 
change the node number of the node itself and any elements that reference it will be updated, but you cannot 
change a node number from the element properties like you can for a primitive surface. 

Element shapes are defined by the nodes. Move a node to a different location and the shape changes. 
Material Orienter may be needed for solid element to define axes for anisotropic material property 

Options exist under Modeling Tools that are useful for working with elements 


Refine Elements: subdivides elements to convert a quad element into 4 quads based on the midpoint of each edge. Similarly this can 
be applied to triangular elements and creates 4 tri elements. 

Reverse Connectivity: At times, the top side of an element is defined incorrectly. The connectivity can be reversed to switch the top and 
bottom sides 

Shift Connectivity: If using edge contactors with elements, it often assigns the first, second, etc. edge for making contact. The marker 
is displayed along the specified edge(s). However, instead of having to guess which edge is the first, the connectivity can be shifted to 
redefine the first edge. The marker will visible shift and this command should be repeated for each element until the markers line up 
along the desired location 

Convert FD to FE: Elements cannot be subdivided. A rectangle defined with edge nodes can be converted to elements (FD/FEM 
Network), but a 3x5 edge node breakdown become 8 individual elements 


Options also exist under FD/FEM Network that are useful for working with elements 


Extrude/Revolve Planar Elements: You can extrude a set of planar element along a line or revolve around an axis to define 3D 
elements. The number of elements along the extrusion can also be specified 

Surface Coat: Solid elements do not radiate and cannot be used with a contactor. Therefore, to use these features, a plate element of 
zero thickness is needed. Surface Coating created zero thickness plate elements on all free faces of a group of solid elements. Free 
faces are defined as faces that are not shared between elements. 


BUILD MODEL: Node Discretization 


Common question: How many nodes should | use? 

The answer depends on the design... 

— Thick, high conductivity material structure...fewer nodes 

— Thin OR low conductivity structure...some nodes 

— Thin AND low conductivity structure...many nodes 

Decide as an analyst your maximum allowable temperature between two nodes. (0.5 K, 2 K, 5 K? and where) 

v If two nodes are very close in temperature, they can be merged/eliminated 

v If two nodes are very far apart in temperature, you may need more in between 

Y Do not force two nodes to be at nearly the same temperature with a very high coupling...merge instead 

If conduction Is... 

Y Mostly 1D, consider only discretizing in one direction (along strut, bar, or heatpipe) 

Y Mostly 2D and planar (e.g. thin panel), you should discretize in two directions 

Y Mostly 2D and non planar (cylinder, cone), consider uneven loading from environment or radiation when 
deciding on circumferential discretization 

Y Mostly 3D, use solids or solid FE 

Often it also comes down to a matter of convenience for modeling conduction across an interface 

After you have run your model for the first time, investigate if the gradients warrant more or fewer nodes 

— Best to stay with surfaces and edge nodes if possible while determining this. 

— It is easy to renodalize a surface; it is far more difficult to remesh elements (unless you use the TD mesher) 

— Can Refine Mesh in TD which makes an FE Tri into 3 Tris and a Quad into 4 Quads...cannot go back though 


BUILD MODEL: Assigning Properties 
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Thin Shell Data 


¢ Once surfaces/solids are geometrically defined, remaining SS 
properties need to be defined in ThermalDesktop cee ; — ; 
VY Subdivision*: X,Y,Z breakdown and subdivision, Edge vs. Centroid — = = Ya 1 n 
Vv Numbering*: Submodel and Node Numbers (multiple submodels allowed) | ~==2 —* OO ——————— 
Y Radiation: Active Sides, ThermoOptical Properties, Overrides =_—- —_= eel wine wee 
Y Active sides defined for each Radiation Group Mest alee: & 
Y Cond/Cap: Thermophysical Material, Thickness = _ 
Y Insulation: Top and Bottom Sides, Insulation Material, Overrides 
Y Contact: Do Not Use eres = 
VY Surface*: Dimensions, Comment = aa 
Y Comment for Finite Elements (not Surface) — =a 
Y Trans/Rot: Additional rotations and translations aan ma 
* These properties not applicable to finite elements ii ea selina tl 
Thin Shell Data Thin Shell Data Thin Shell Data 
pret eer gy gg corer a 


Use List: 


Free Molecular Conduction Accommodation Coef from Emissivity 
Top/Out DEFAULT 


Bottom/In DEFAULT 


Note 1: Insulati 
Opti 


EFAULT 


Top/Out Side Node Numbering/Creation 
Offset Node ID's b 10000 


1 Thickness: 0 


Muttiple Materials (Stack) 


Bottom/In Side Node Numbering/Creation 
© Offset Node ID's by: 200000 


Use new submodel 


Calc Type: [Based on material property 


Init Temp: 20 c 


Overides... 


.a changes due to changing radius are not accounted for. All inear calculations are K"A/L 


jon conductors are assumed constant area. Arei 
Note 2: Optical Properties of the Insulation are set on the Radiation Tab 


BUILD MODEL: Node Numbering 13 


* Each node is an independent calculation point in SINDA [ Merge Coincident Nodes 
Y Use duplicate numbers if you intend for the temperatures to be the same (e.g. @ I/F) 


* If node numbers are duplicated and nearly coincident, best to merge them into single node entity decide a 

¢ Use submodels for organization to identify component associations (e.g. Opt_Bench, Scan_Mir, FPA) ace eer [Fist Selected 
Y Too many submodels gets cumbersome in the ModelBrowser...don’t go overboard Smallest Node 

¢ But what node numbers should be used? ~ Lesser Submodel Name 
— Good practice to number nodes according to their location within an assembly : 
Y Best to keep MLI offset at 100000 for one side and 100000/200000 for both $$ 


Y Number first portion of a component in the 1xxx range, the next in 2xxx, and so on... = = a = 
— Could start with 1xxxx range, but run out of “blocks” sooner before the 100000 range [sou | 1 eee ll 
— If you have more than 1000 nodes in a component either skip ranges (e.g. 1000-2999, 3000-3999, 1° at Saco 
4000-4999, etc) or use a block of 10000 for the entire component i im | of permite 
¢ Why bother? Why not just renumber everything and let TD display the contours? sae | aca 
— Helpful for navigating model to have ranges...If 1000-1999 is the motor shaft, can turn its visibility I a i Pater 
on and off in the ModelBrowser by selecting that range + ine | 4 ocala 
— Helpful for post processing...if all 100000 nodes are MLI, can likely ignore their temperatures. If eo | are 
1000-1999 is the motor shaft, can apply limits to those nodes knowing what they represent tow | pil Fenstmuct 
— Helpful for documentation...likely you will not be the last user of a model. The next person who a | 
picks it up and has to work with it will have an easier time learning the model if it is organized Por 4 peo 
* To propagate numbering scheme in Z direction... z amt : PT ee 
v Renumber all nodes to 100000 a: ea) 28>: paar 
v¥ Renumber first layer of nodes to 1000 reer. =I Ee Q 2 ! : foro 
Y Copy 1° layer of nodes in z over top of 2" layer Seis refuse 197 - om} ree 
v Merge nodes and keep Smallest Node ID fete ante eoerere Simms: Cains 
v Renumber with Add Increment checked foam! ay TEerusesreh 
* coi 2019 © <| | | cs) TELE FOAS HEATER 
Y Repeat [tance] (Help a0" sbjects selected Jz | Hotning selected 


BUILD MODEL: Radiation Groups 


¢ While Thermo-optical properties are specified for a surface, it does not necessarily 


mean that the surface participates in the radiative calculations 
¢ Radiation Groups define collections of surfaces (i.e. enclosures). Only one 
Radiation Group is considered at a time during calculations 


¢ For each Radiation Group, participation in calculations (activity) can be specified 


¢ Multiple Radiation Groups can be used to define different configurations 

¢ Radiation groups can be OR merged (If Top is active in on and Bottom is active 
in the other, OR merging makes Both sides active) 

¢ They can also be copied to spawn a different configuration which may be used 
for changes in the geometry such as use of different instrument models or 
removal of a cover 

Y Define the Radiation Groups before assigning activity 

Y For each surface, assign whether the top side, bottom side, both sides, or 
neither side participate in the solution. Additionally, Not in Group can be used 
for surfaces that should be excluded completely from calculations 

¢ Active None vs. Not in Group: Active none will still reflect and absorb energy; 
it just won't be tallied for the surface. Not in Group: surface is not even there 


Y Once activity is assigned, visualize using Model Checks-Display Active Sides. 
This will be for a specified Radiation Group. Colors indicate activity: Yellow ( 
(Both Sides), Green (Visible Side), Light Blue (Not Visible Side), Dark Blue 
(Neither Side), Red (Not in Group) 

Y Be sure to check activity for every group that will be used in calculations... 


Radiation Analysis Group Manager 
‘CORONA Add 
INT_AMS a 
INT_IC 
INT_PM Rename | Purge Unused | 
INT_SM 
INT_\WFI_IFM opy 
WEFIRST_External [default] 
WFIRST_External_Orb 
acte 
Bien 
Thin Shell Data se | 
Radiation | Cond/Cap | Insulation | Comment 
Analysis Group Name, Active Side Optical Properties for Radiation Calculations 
[CORONA vg} 
INT_AMS ry 5] Top/Out WFI_GBK = 
INT_IC [n/a] 
INT_PM In/a] Bottom/In DEFAULT ba 
INT_SM [n/a] 
INT_WFI_IFM [n/a] 
WFIRST_Extemal top/out = 
WEIRST-Edemel Orb top/out | Top Side Overides.. Bottom Side Overrides... | 
Edit Active Side 
Top/Out lommodation Coef from Emissivity 
Bottoméln oa 
Both oa 
None [will reflect and absorb energy) 
® Not in analysis group (Not part of the calculations) 
[Cancel | 
Edit... 


is group: | WFIRST_Extemal | 


BUILD MODEL: Controlling Visibility 1s 


Seldom is it useful to display the entire model. Most often, only one part of the model is being updated at a time 

Visibility control allow only parts of the model to be show. Three ways to control visibility: Layers, Global Entity Visibility, ele Level Visibility 
Layers: (layers command, native to AutoCAD) re ee 
— Ifa layer is Off or Frozen (Blue Light Bulb or Snowflake) anything on that layer will not be displayed sock a en ent init Site 
— Also can control from drop down along toolbar 

— HINT: If something isn’t showing up in post processing, make sure the layer is not frozen in the VP too! 
Global Entity: (Thermal-Preferences) 


— If the visibility of a particular object type is unchecked, it will not be displayed 


oo BUS_Avionics 

of Ol BUS_Deck-Bot 
of Ol BUS_Deck-Top 
of Hl BUS_HeatPipes 
of BUS_HGAS 

of O BUS_MidStructure 
of Dl BUS_Module7 
of Ml BUS_Prop 

of Dl BUS_Radiators 
of Ml BUS_SA 
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# 
RB 
4 
iu 
pitty 
a 
ai 


Patt d 
OF: be 
roll » 
m $5 
2 38 
° ces 
9 
Y 
Ox 
P - @ kirt 
— On/Off Toggles for common objects (Node, UserNode, Surface, Solid Element, Heat Load, Conductor, © 38 aim ars | 
. Dy , Sunshield 
and Contactor) available on toolbar rE 9 8 em N00 spon 
. 7 G CGLBENCH 
Object Level: (Blue and Yellow Lightbulb toolbar buttons) = x os eae ca omen 
- . Peer wen P 9 2) of CGLDM1 
— Objects can be selected in Model Browser and visibility turned on or off (Blue and Yellow Lightbulb) ale 9m CoLDMe 
F Step is exad L 
— Objects can be selected from the screen and visibility turned off |e gecdmcarae 
— Objects can be selected by group (even if currently off) and turned on alc 38 ato canine ca nous 
— On/Off state can apply to solid CAD but hard to manage; use “all” for selection if needed ai" eeu cecum 
Best to use Layers for CAD objects. User choice for Thermal objects, as visibility can be controlled by Model Browser 
a] niin iaiaiaininiemeaiuaaiiaies fesse = = Visibilty | Graphics Size | Curved Elements | SINDA | Advanced r ma wi a wv EXkv ae 2 - 
List Edit Display Options Help * & Filters «|[s. Name ae Trans Plot S.. Pau Nu Description “ 
i] [IR & KOK AL) [A EVEe, Ie (>) (x) bg I }TD/AC Nodes! ¥| Lumps | Conductors Nodal aay : °S ae _ regen ; - : - . = . 
la BUS_Deck-Bot & 0 lr2 BR 
Submodel Node Tree a (¥] User Defined Nodes | Paths 7| Contactors 7] Edges in Shaded Views > BUS Deck-Top ; a 0 jor? & 
; Bison li owes ike ewes |) 7 ca 6B 
=) 1 [9 Solid Finite Elements [¥] Pipes ¥] Material Orienters 4 BUS Midstructure Q 0 lo2 BB 
=} BOUND_MAIN.1;:12EEE = Wi ibcseack Se ee ws cas ler os g i : | : a4 
OT Cylinder-Blocker Surface[BOUND_MAIN}::12EED mines ~ 7 pabaatieely . x + ace 
z BOUND_THRUST ) Mesh Sapukiis ¥) FTies ¥| Assemblies Pipe Surfaces Drawn |= BUSSA Q ts 0 lors @B 
:3-(H] BUS_BATT Z)iFoces peel :* ) jonae 
+ {st] BUS_CDH v) ee Cain {¥] Heat Exchangers ¥] Contact Conductance a BUS Sunshield g om 0 I eg 
fl BUS DFCE coe M — oem: lee . fee 
+ BUS_DIPLEX aii ae aa 3 ps CGLBENCH is Q # 0 | 8 8 
+)-{sM] BUS_EPC — — | Color Contours 2 l= CGILCMR_ELEX 9 & 0 lors BB 
21 fi] BUS_HGAS — — i Pn ee . wee 
+)-{si] BUS_K_BAND_MOD 3 | Chinen tine all ages ot ar 
eH BUS_MACE A All: i layers displayed of 118 total layers 


BUILD MODEL: Defining Variable Parameters ° 
¢ Symbols can be used to control just about anything in Sprbl anaes 
ThermalDesktop 7 a 
— Dimensions (thickness, length, width, height) ———— = ——— (ia 
— State Flags (On/Off, Enabled/Disabled, Open/Close, Stow/Deploy) Meni es shikai : 
— Margin Factors (Heat Load Multiplier) [rein 
— Other Values (Conductance, Heat Load, Setpoint, Rotations, etc) tes See epee more 
¢ Symbols can | 
— be single values (Integer or Real), expressions, or arrays of values oe — 
— depend on other symbols Spacecat Iteiace Temperate Alas Output Thermal Beskiog Symbal as SGA Regt 
— be overridden in the Case Set Manager for a particular analysis . ee 
Ps Symb ols should wiaimanies = @ cutout Symbol Expression 
Vv Be named intuitively. Use prefixes for thickness (thk_), heat (Q_), |__, — smsmsowsesnoanens. acre 7 
conductance (G_), Multipliers (Mult_), etc. Fis der usp dra mv gota Gots) (econ 
Y Be grouped together where logical (power, subsystem, etc) Ee 
Y Include a description or comment aa = 
* Symbols can be passed through to SINDA as REGISTERS ES lai 
— This is automatically done if the symbolisneededinSINDA = = |. a . 
— User can override to ensure symbol is passed as REGISTER if not 
detected by ThermalDesktop ae 
¢ To use symbol, dbl click in textbox to access Expression Editor . 
* Values can be set in Symbol Manager, Case Set Overrides, or ee ee 
user defined logic...be careful where you define things a Se 


BUILD MODEL: Assigning Variable Parameters 


¢ Once a symbol is defined in ThermalDesktop, double click in just 

about any textbox to access Expression Editor to use it 

— Values that are dependent on expressions are indicated in Bold font 

— Values with a description are highlighted in light blue 

— Illegal or incalculable expressions default to 1 with a red background 
(importing a model without importing symbols) 

Unit conversion can be performed 

Can output expression to SINDA (as REGISTERs) 

Conditional statements (Not Valid in SINDA FORTRAN block...) 

— (Heater_ON ==1)?10:0 

— If Heater_ON is equal to 1, set expression to 10 otherwise 0 

interp function can be used to access elements of symbol arrays 

— interp(CaseArray,PowerArray,CaseFlag) 

— Look in CaseArray for Case Flag and return entry at same index in 
PowerArray — Not Valid in SINDA FORTRAN block... 

Example symbol usages include: 

— Power Dissipation (Heat Load or Heater) 

— Conductance Values or Convection Coefficient 

— PID Gains 

— Assembly Rotation Angles or Translations (Open/Close) 

— Thickness 

— Geometric Distances 


Expression Editor 


Select units for: Temperature 
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Expression: 
Temp_SC_IF 


Add Symbol... 


Node - Multi Edit Mode 


Enabled... 


Submodel,  VIS_DET_MTG 


Comment: tial temp: 293 K 
Set Spacecraft Interface Temperature 
¥ Output Above Expression To SINDA (Please make sure that all REAL values have a‘. with the number) 
Expression is in SINDA Units (No units Conversion performed) 
Disable Warnings for this Expression 
|__ Cancel Help 


Boundary Temperatures 
Initial Temperatures 


— Setpoints 

— MLI Assignment 
Number of Rays or Error for MCRT 
Margin Factors/Multipliers 

State Flag (On/Off) 

Failure Flag (zero out value to disconnect or null) 
Absorptivity/Emissivity/k/e*/etc 
= Fic. 


BUILD MODEL: Moving Geometry 


* Two methods exist in ThermalDesktop for handling moving Edt Assembly Edt Aasembly = 
ke Assembly | Trans/Rot Assembly | Trans/Rot 
geometry: Assemblies and Trackers 
Name: Linnamed| Translation X: in 
— For both, geometry must be attached nen oes . ; 
— lf using FE, nodes must be attached since nodes define elements ~- . oe ; ’ 
e@ AS semblie S Graphically Display Name Rotation 1: 0 |X y| Degrees 
_ Everything attached moves when assembly is moved ¥| Display displacement vector and base coordinate system Rotation 2: 0 |Y y| Degrees 
‘ é é. Rotation 3: 0 |Z | Degrees 
— Be careful using AutoCAD move command; moving TD objects and the Y) Active — 
assembly makes everything attached to the assembly move twice 
— Useful for moving geometry based on user needs (e.g. door opening, shield — Coc} [_cree_} Le] 
deploying, lid closing, etc). Could be orbit position dependent. pe arererrereane eae = 
e@ Tracke rs 7) Do Fast Spin Calculations Name: unnamed 
Symbol Name: | TEEEWINE! M Track Working Mode 
— Similar to Assemblies, but point specified axis towards target by rotating © tes Ean tomes @ Sun © Amys 
about second specified axis. Can be nested, but gets complex quickly ane = in eed 
Number of Increments: 12 sal fk =A 
— Useful for simple solar arrays or more complex spacecraft pointing se eget Fight Ascensior: (0 Beret eer (teal 
— Utilized by HR and Articulating Radk calculations UO | aaa S z 
¢ Fast Spin (many times in one orbit position) tre Tue Anema 
. A op. < 7 % . Program Mean Anomaly 
— Fast spinning capability exists for radiation calculations for scan mirrors, f= as Disable (Never Working 
. ® || CAESAR_GCS_SCS_Oct6A_PCM_Final5.dwg 
rotating drums, etc. List Edit Display Options Help Program _ | 
— Use an assembly, assign spinning geometry, and assign rotation angle to my eee es Dis 
sym bol Va riable o 1 hae —— ; Size: 1 in 
: : ; ae 5. J. Articulator-GCS_Hinge::1D0C | To: 360 dearees 
— In Radk or HR Spin Tab, specify start angle, end angle and increments ge 
H _— ¥) Active 
— ity, -- ). Group 2 
Be sure to use 0 error and max number of rays per position : oe Sees se ie 
¢ Slow Spin/Slew (<< one time in orbit position) Fl Gyindef GCS F PLATER20616 
a . a © GCS_IF_PLATE.1001::2971B 
— Orbit Rotation or slew to multiple of hrMeanAnom © GCS JF PLATE1002:29717 


BUILD MODEL: Establishing Connectivity . 


¢ When nodes are not merged to establish connectivity between two surfaces, a user must add this manually. Two 
methods exist in ThermalDesktop for this: Conductors and Contactors 


Conductor: Tie a specified node to a 
surface, multiple nodes, or a single node 


i Conductor 
| Enabled... | 
Comment: Hinge between GCS and GCS IF Plate 
Submodel: GCS_STRUCT a 


Best to use for 


® Auto-number ID 


abel Point-to-Point or 
Type: Genetic H p. s | 
ee to-Surface 
Area/length: 0.055 mm ofitact 
Vs. Temp Diff 
¥| Use material: y 


Radiation conductor 


Per Area Bi-directional Heat Flow M4 | 


Insulation Nodes 


From Node: GCS_STRUCT.23061::27F6D | Reselect | 

To (Uses Area): 

Quad Elem[GCS_IF_PLATE]:297D2 Top 

\Quad Elem{GCS_IF_PLATE]:297D3 Top fe se ] 
al. 8 |O/A 
tvs E@, 


Cancel | Help 


¢ Both 


— Can be material dependent, time or 
temperature dependent or constant value 

— Can be to insulation nodes, one way, or 
radiative conductors 


¢ Conductors 
— Remeshing/renodalization can break this link 
— Can be per Area, Cannot mix types (i.e. cannot 
add a surface to a node to node conductor) 


* Contactors 

— Renodalization does not break this link, unless 
new elements are created 

vY Must specify Edge or Face contactor type and 
edges/face sides to make contact 

Vv Use Ray Trace for Face Contact 

Y Specify reasonable tolerance 

v Preference is to use Per Area or Length over 
absolute conductance 

Y Often, best not to Apply Surface Thickness 

— Contact From displayed as Green Arrow, To 
faces as Yellow Arrows 

— Can visually display contact made (or missed) 


Contactor: Tie nodes associated with two 
surfaces or two edges together 


r 
Contactor (=| 
Enabled... Add Code For Conductors 
Comment: Bottom Tank PCM Core to Housing a 
Conductor Submodel GCS_PCM +’ 
Contact From: Faces Y | Restart Files From Current Default Directory Y 
Conduction Coefficient (Conductance/Area] 
0.00155 Wémm*2/C [_ Amray... Vs. Time Vs. Temp Diff 
Input Value Type: | Per Area Or Length ’ 
Use Materiak =| DEFAULT 
raisin lew BeSt to use for 
Use U Scaling Ed e O r 
Use ¥ Scaling v T 
Use Insulation Nodes Su rfa ce-to-Su rfa ce 
Inputs For Connection Algorithm 
Integrations Intervals: 10 Point onta ct 
Tokens 5 a Max Check Objects: 
Apply Surface Thickness To Test Points © Ray Trace Algorithm 
Oct Cell Subdivisions: 6 
From (5 objects, Area = 63039.93 mm2): To (5 objects, Area = 69725.6 mm2): 
Solid CyI[GCS_PCM]::287F7 HMIN Disk[GCS_PCM_HOUSING]::288B7 
Solid Cyl[GCS_PCM]:28796 HMIN Disk[GCS_PCM_HOUSING]:28775 
Solid CyI[GCS_PCM]:26D28 HMIN Disk{GCS_PCM_HOUSING]::287F5 
Solid CyI[GCS_PCM]::28858 HMIN Disk[GCS_PCM_HOUSING]::28856 
Solid CyI[GCS_PCM]:24172 HMIN Disk{GCS_PCM_HOUSING]::24202 
H(ZIZ Rl (RIRIR Olan. Ee, jo [leh fale ie lola fe. Ee, 
| OK | | Cancel | | Help | Show Cales J 


BUILD MODEL: Establishing Boundary Conditions 


* Boundary conditions must be applied to models. Typical — . — — | 
boundary conditions for thermal analysis include: fixed |“ =" Ane ETE 
temperatures and fixed heat loads came ————— 
— Boundary nodes defined by overriding the calculationsby "| page 

elements/surfaces. Need to also specify temperature a <r 

* Heat Loads eis eee 

V User specified application of heat to nodes, surfaces, a Cine SRSE_S5M) 258 Ou <Ch 
or solids as —— GR I@IO [A 

Y Can be total load or flux eee |e Pa, 

— Can be temperature dependent or time dependent ne LgeSunoiet SDE. BP 

V Specify Value and assign appropriate submodel = hell atad cee 

¢ Heaters Heater Power (5 w Set To Mid Pont Temperature 
v User specifies application of heat location and sensing location tee a Sone 
— Can be applied to nodes, surfaces, or solids OfTenn: 220 K ce * 

V Specify Available Power, On Temperature, and Off Temperature a ee 
Y Specify control approach (thermostatic or proportional) ( Sereetes._| (Foon) (Bese Meta. 


Y Specify Steady State behavior (Proportional recommended) epee 
Y Specify Append string for symbol output recognition Rae DET WTGAIOTTSOSS Top 
¢ TD creates power, total power, on time and #cycle registers 
using this string 
— Can specify sensing method (defaults to area weighted 
average temperature, but could be min temp or user defined) 


S rature From: 
VIS_DET_MTG.1017::1B67 Top 


lhe fe) RIL JO lA Be Ee, le 2 Te Fe] RIO |OJA Bee, 
| Cancel Help ] 


BUILD MODEL: Modifying Entities in Objects and Domain Tag Sets 


21 


The entities referenced by Heat Loads, Heaters, Conductors, and Contactors can be modified from those references when 


originally created 
e Add (+), Remove, Remove Selected, Edit 


ele le fe) ILO lA be Ee, 


¢ Switch Visibility, Turn on Visibility, Turn Off Visibility, Show Only Selected (O) or All (A), Turn on, Turn off Node Numbers 


Domain Tag Sets are groups of objects that can be 

used in numerous ways 

Create the Domain Tag Set 

Y Specify the Name 

Y Specify the type of object it will contain (cannot 
change this later!!) 

Y¥ Add/remove objects using buttons described above 

Domain Tag Sets show up in Model Browser and 

are useful for controlling component On/Off 

visibility 

Used as Groups for mapping temperatures to FEM 

Domain Tag Sets may be used as target objects for 

Conductors, Heat Loads, Contactors even if it 

does not contain any objects yet 

Useful for defining interfaces that may change on 

either side (dummy plate vs. SC deck) 

When prompted to select objects, type d to bring 

up list box. Only Domain Tag Sets of a valid type 

will be shown 


= 


List Edit Display Options Help 


Pe [F MEI 1@. [0 [A EV.E@ Ib (> I[<) 


1] WFIRST_PreSRR_Discipline Analyses v3M_Thermal_Cl.. [oco-|[-(=)-|[mbsau} | DO™#I7 T29 Set Manager 


Domain Tag Set Objects 
2 BUS_MAP_BOT_DECK 
3 BUS_MAP_CENT_CYL 
9 BUS_MAP_CORNER_LONGERONS 
9 BUS_MAP_HGAS_BOOM 
{2} BUS_MAP_HGAS_BRACKET 
2} BUS_MAP_HGAS_DECK 
2} BUS_MAP_HGAS _DISH 
9-3 BUS_MAP_HGAS GIMBAL 
4) BUS_MAP_HGAS HINGE 
{23 BUS_MAP_HGAS_LOCKS 
5 BUS_MAP_IC 
-{ BUS_MAP_IC_MONOPODS 
2 BUS_MAP_MODULE1 
2 BUS_MAP_MODULE2 
2 BUS_MAP_MODULE3 
2} BUS_MAP_MODULE4 
4 BUS_MAP_MODULES 
4 BUS_MAP_MODULE6 
4 BUS_MAP_OBA 
{2 BUS_MAP_OBA_BASE 
{2 BUS_MAP_OBA_FCR_SUP 

GS 


1 BUS_MAP_PAF 

3 BUS_MAP_PROP_DECK 

9 BUS_MAP_PROP_STIFFENER 
2} BUS_MAP_PROP_STRUCT 
2 BUS_MAP_RADIAL_PANELS 
ce)-{23 BUS_MAP_RWAS 

2} BUS_MAP_SA_MID_MX 


£2 mr A RAT mw 


CanceledFiltering from 5 objects 
hw objects meet criteria 


Select surfaces to remove or [Domain_tagset/MB/GRP]:*Cance 


" 


Domain Tag Set Name Object Type 
’ BUS_MAP_BOT_DECK [ANY SET 
BUS_MAP_CENT_CYL [ANY SET Create 
BUS_MA4P_CORNER_LONGERONS [ANY SET 
BUS_MAP_HGAS_BOOM ANY SET Fienome 
BUS_MAP_HGAS_BRACKET [ANY SET C 
BUS_MAP_HGAS_DECK ANY SET apy 
BUS_MAP_HGAS_DISH [ANY SET 
BUS_MAP_HGAS_GIMBAL [ANY SET 
BUS_MAP_HGAS_HINGE [ANY SET 
BUS_MAP_HGAS_LOCKS [ANY SET Delete 
RIC MAP IC [ANY CETI 
Contents of the Selected Domain Tag Set: 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABCA a 
Brick Elem-Honeycomb Core[BUS_STRUCT)::ABCS + |Z 7 ] 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABD2 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABC4 ny Is i io iA 
Brick Elem-Honeycomb Core[BUS_STRUCT)::ABC2 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABCO * 2%, 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABCF 
Brick Elem-Honeycomb Core[BUS_STRUCT)::ABCE 
Brick Elem-Honeycomb Core[BUS_STRUCT}:ABCC 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABC? rn 
Brick Elem-Honeycomb Core[BUS_STRUCT]:ABBC Select Domain Tag Set 
Brick Elem-Honeycomb Core[BUS_STRUCT)::ABBB REE — = | 
= = —|| BUS_MAP_SA_MID_MX 
CGI_BENCH 
CGI_BIPODS 


Domain Tag Set Name: 


Objects Type: 


2 Set 
{ OK Node Set 
Solid Set 


>.|~Select surfaces to remove or [Domain _tagset MB GRP]: Any Set 


te] I | [|] 2 Jk] =| + | 


Create New Domain Tag Set 


DOMAIN_TAG_SET_1 


CGI_FRAME 
CGI_HIFED_BENCH 
CGI_MAP_CAMERAS 
CGI_RADS 
NOBA_MEASURES 
SC_TOP_TECK 
TOMA_BENCH 
TOMA_BIPODS 
TOMA _F1 
TOMA_F2 


TOMA_M3 
TOMA_M4 


BUILD MODEL: Simulating Thermal Hardware 


¢ Every analyst will develop their own modeling style (i.e. how physical components will be represented in the 
analysis model) 
¢ The following slides include suggestions on methods to model: 


¢ Honeycomb Panels 

¢ Insulation 

¢ Thermistors and Thermocouples 

¢ Constant Conductance Heat Pipe 

¢ Variable Conductance Heat Pipe 

¢ ThermoElectric or CryoCooler 

¢ Controller (most often for a heater) 
e Phase Change Material 


e Louvers 


BUILD MODEL: Honeycomb Panel Modeling 7 


¢ Three ways exist to model Honeycomb Panels 

1. Neglect Core: use a single sided surface with 2x the 
facesheet thickness and the material of the facesheet Mutipiers: 

2. Include Core through Conductivity: use a double sided Density: = 1 
surface with thickness and material specified for facesheet Yor — ' 
core, facesheet. Note: the in-plane spreading effect of the "Yo" | 
core material is not included using this method... — © 

3. Include full core effect: model the core as a solid element 
and the facesheets as plate elements using the same nodes 
as the top/bottom for each facesheet. If applying non- 
isotropic properties, must define material orienter for solids 

¢ TIP: you can use the U Cond/V Cond multipliers increase the 


conductivity as a function of the core material to account for 
inplane spreading by the core —— 


Comment: Core Hements for Honeycomb Panel 


HC Core Material Orienter ~~ AL_S06ITE. 3p Cr 
(defines CS for anisotropic materials) | Mssisoiet— He cae 


Cond submodel: PANELS 


Multiplication Factors: 


Density: 


1 

» Conductivity: 1 M et h O d 3 
1 
1 


Y Conductivity: 


2 Conductivity: 


Cancel Help 


eT 


Subdivision | Numbering | Radiation | Cond/Cap | Contact | Insulation | Surface | Trans/Rot | 


[___ Generate Cond/Cap | Method 1 

Cond Submodel: PANELS ’ 

Gen Nodes: —_| Based on material property y| 
Material Thickness(n) 
M55J_Composite > 0.08 


DEFAULT 0 


DEFAULT 0.0393701 


Ragiation | Cond/Cap | Contact | Insulation | Surface | Trens/Rot 


-_ Method 2 


Submodel: 
PANELS 


@ Use Stat ID: 1 @ Use Start ID: 101] 


© Use List: Use List: 


Material MI et h O d Z Thickness(in) 
Top/Out: M55J_Composite ¥ 0.040 
Separation: AL_6061T6_3.1pcf_Core ~ 0.500 
Bottom/In: M55J_Composite > 0.040 


Perey 


BUILD MODEL: Insulation Modeling 24 


¢ Insulation may be conductive or MLI (i.e. radiative) 


Y Define insulation material 


Vv If material has conductivity <> 0, k and thickness will be used 
v If material has s* defined, radiative will be used (thickness ignored) 
Y Assigning kg/m3 of 0.6 and applying a thickness of 1 m will allow MLI 


mass to be estimated by model 
— Possible to have both... 
¢ Use intuitive node numbering 
Y Keep all non MLI nodes below 100000 


Y Use Node Offset of 100000 (200000 if both sides have insulation) 
— When post processing, any nodes above 100000 are MLI temps and 


can likely be disregarded 


underlying surface 


the P button at the top 


isolation to have low ¢€ and low ¢*% in series 


Thin Shell Data 
Subdivision | Numbering | Radiation | Cond/Cap | Contact | Insulation | Surface | Trans/Rot 
Put on top/out side! (FJ Stack Manager | [WJ Put on bottom/in side [2] 
Top/Out Side Material/Thickness Bottom/In Side Material 


® Single Material 
Material: = e"nom 
Thickness: 0 m Number of Nodes 1 
Multiple Materials (Stack) 
DEFAULT 


Top/Out Side Node Numbering/Creation Bottom/In Side Node Numbering/Creation 
Offset Node ID's by 1 @) Offset Node ID's by: 200000 


Use new submodel. 


Based on material property Calc Type: | Based on material property y| 
Init Temp: 20 c 
rs are assumed constant area. Area changes due to changitt="=s @Paccounted for. All linear calculations are K"A/L. 


D eterPConduct 
Note 2: Optical Properties of the Insulation are set on the Radiation Tab 


If overrides are used (i.e. applying MLI to only some of the 
nodes), make sure to adjust the optical property overrides as well 
Y Opt Prop overrides will need to reference MLI node number, not 


MILI can be “programmed” to either be enabled or disabled using 


For “tented” MLI (offset by large gap from surface), best not to 
use low € optical property for inner or outer layer — doubles 


Thin Shell Data al 


ey) 
Subdivision | Numbering | Radiation | Cond/Cap | Contact | Insulation | Surface | Trans/Rot 
Analysis Group Name, Active Side Optical Properties for Radiation Calculations 
|EESS both 
ane bath | Top/Out: Z93P_WhitePaint - 
Bottom/In M55J ¥ 


Free Molecular Conduction Accommodation Coef from Emissivity 
Top/Out: DEFAULT oa 


Bottom/In: DEFAULT oa 


Edit... 


BUILD MODEL: Thermocouple and Thermistor Modeling 


¢ Thermocouples, Thermistors, and other measurement devices are the tangible data 
readings visible to the thermal engineer based on the actual hardware 
¢ ThermalDesktop has an object type called a Measure, which represents a spatial 


location for a temperature reading. 


¢ Measures are displayed graphically with their Name and show the mapped location. 
They can be moved to the exact location desired by the user 

Y Define the Measure and place it correctly in space 

Y Set the size large enough to be visible when displayed 

Y Assign a meaningful name (e.g. Telemetry mnemonic) 

Y Determine if you wish to output measure value to an existing 


Register or a Boundary Node 


Y Good idea to Create groups to ensure Measure maps to 


reasonable entity 

V Set tolerance 

v Execute mapping 

v Snap Measures to 
Mapped Entity 

— Can copy measures from 
One drawing to another 

— Can import from text file 


is 
Variable Tolerance 


Input one value per line, in units of mm 


Ignore Thickness when Gen Cond/Cap is not enabled 


‘Use Advanced Mapping 


Cancel 


N 
uw 


Edit Temperature Measure 


Sinda Interface 
Output Register 
Register Name: 
¥| Output Node 
Submodel: NOBA_HTR 
ID: 2599 


Use Conductor and Thermal Capacitance 
Conductor Value: 1 wik 
Thermal Capacitance: 1 IK 
Connection to Model 
Test All TD Entities 


@) Test AutoCAD Group: NOBA_MEASURES 


Connect to Outermost Nodes of Insulation (if found on surface) 


Lx) |_Cancel _] 


BUILD MODEL: Constant Conductance Heatpipe Modeling 


¢ Heatpipes are commonly used devices in thermal designs 


— They operate on an evaporation/condensation cycle that allows heat to be 
transported over large distances with a minimal temperature gradient 

— Heat enters the evaporator end of the pipe and vaporizes the working fluid 

— The vapor is driven to the condenser end of the pipe by a small pressure gradient 

— The vapor changes to liquid in the condenser end releasing the heat 

— The liquid is pumped to the evaporator by capillary forces through a wick structure 

— Unfortunately, the capillary forces that drive the fluid are very small compared to 
ee which places constraints on HPs during ground testing 
Device can operate in reverse depending on loads... 

HEATPIPE exists in FloCAD, but a simpler way exists to model a heatpipe that 

does not require a FIOCAD license 

Y Define surfaces based on length/shape of pipe 

Vv Set as edge nodes 1 x Z breakdown, assign nodes to increment from one end to the 
other (e.g. x001, x002, x003...) and merge at interfaces 


Y Set width of surface to width of interface flange 

Y Set Material to Aluminum (Assuming aluminum walls) 

Y Set Thickness to effective thickness (HP Cross Sectional Area / Width) 

Y Create Arithmetic node to represent vapor node and set to (x999) 

VY Create a Node-to-Surface conductor from Vapor node to all HP surfaces 

Y Check the Per Area option and set conductance value to h per linear inch / Flange 
Width (3.0 W/lin in. K typical for ammonia HP) 

v 


Add contactor from evaporator/condenser sections to structure/radiator 


Node 


Enabled... 


Submodel: HP 


ID: 1999 


Comment: Heat Pipe Vapor Node 


Initial temp: 20 iE 


Type 


Diffusion 


DEFAULT 


®@ Arithmetic 


Conductor 
Comment: Heat Pipe 
Submodel: HP 


®@ Auto-number ID 


ID number: 0 
Type: | Generic 


Value: 1860 


Use material: bolt 
Radiation conductor 
¥) Per Area 


On Nodes 


Dim? 20 z= 


| Bi-directional Heat Flov 


From Node: HP.1999-Heat Pipe Vapor Node::AF2 


To (Uses Area): 

Rect-Heat Pipes[HP]::AE7 Top 
Rect-Heat Pipes[HP]::10FB Top 
Rect-Heat Pipes[HP]::10E4 Top 


BUILD MODEL: Variable Conductance Heatpipe Modeling 7 


¢ A Variable Conductance Heat Pipe (VCHP) is an extension of a CCHP 

— Itincludes a reservoir at the condenser end filled with Non-Condensable Gas (NCG) 

— This reservoirs is cold biased and includes a heater to maintain its temperature 

— When the reservoir heats, the NCG expands and partially fills the condenser section, blocking 
condensation and reducing the effective condenser length 

— When the reservoir cools, the NCG contracts and recedes toward the cold reservoir, increasing 
the effective condenser length 

— By heating and cooling of the reservoir, the evaporator end of the VCHP can be controlled to a 
given temperature under a variety of source and sink values 

HEATPIPE routine in SINDA can accommodate modeling a VCHP, but requires additional 

inputs: [Vapor Node, Array of Vap-to-Wall Cond Numbers, Array of Segment Lengths, Volume, HP 

Diameter, Foongs Neyap» Res Volume, Res Temp, NCG charge, Working Fluid and NCG IDs] 

Y Define a heatpipe as previously specified and ensure that the nodes increase incrementally from 
the Evap to the Condensor. Ensure the vapor node is arithmetic! 

Y When defining the Vapor to Wall Conductors, be sure to specify the ID and not Auto-number 

Y This ID will be a fixed value set by the user to be specified in the Vapor-to-Wall Cond Array 

Y Besure to add a geometric representation of the reservoir and its radiator and heater: heat 
should be applied at the Res Radiator and sensing at the evaporator 


Conductor 


Y Define Arrays with Segment Lengths and Cond Nos in the Logic Manager 
Y Define a VARIABLES 1 block in logic manager with a call to HEATPIPE and add in | cm: ——_ Heat Pree 
the appropriate variable and argument values Submodek HP 


® Auto-number ID 


VY Refer to SINDA/FLUINT manual for specifics...be careful with initial conditions ibinaneer 


BUILD MODEL: TEC or Cryocooler Modeling - 


¢ ThermoElectric Coolers utilize the Peltier principal to induce a temperature gradient between two junctions 
when a current is applied to provide cooling. Typically about 6% efficient 
Performance based on 4 related variables specific to device. Knowing 3 can allow 4" to be calculated 
Hot Side Temperature (can be retrieved from model) 
Cold Side Temperature (usually the goal or setpoint) 
Cold Side Load (can be calculated from model) 
Hot Side Load (including input power, current, or voltage), often the independent control variable 
User must decide how to model controller... 
— Assume it can achieve control temperature and set cooling point to boundary temperature 
— OR Model controller with feedback and apply negative cooling load (see PID controller slide) 
If assuming controller can achieve temperature... 
Y Need to extract heat removed from boundary cold side (QFLOW, QFLOWSET or HNQCAL) 
VY Retrieve Hot Side Temperature from model 
V Determine cooling DT from T,,; - Testa 
V Look up power needed for TEC to remove Q heat, and achieve DT for current T,,; 
Y Apply power to Hot Side node 
Y Routine developed for WFC3 to characterize TEC performance curves by 4* order polynomial 
Y Could also use TRIVARIATE array with enough data points... 
¢ If constant power TEC... 
V Determine heat removed based on input power, DT, and T, 
Y Apply a negative heat load at cold side node 
¢ Hot Side and Cold side usually not coupled unless TEC is off... 


BUILD MODEL: Controller (Heater or Cooler) Modeling 


¢ More thermal models are asking for controller modeling: input 
value is mathematically evaluated to predicts output value of 
controller (e.g. Input: Temperature, Output: Heater Power) 
* Common example is thermostatic or proportional heater: 
Y Define On and Off Temperatures and Max power to apply 
Y Define Location to Apply Heat and Location to Sense Temperature 
Y Use Proportional for Steady State behavior 
Y Determine if Proportional is needed for transient 
Y Often needed for tighter stability 
VY Widen Range between Toff and Ton for better control 
¢ For a full PID controller, use a PID controller under Logic Manager 
Y Define Proportional, Integral, and Derivative Gains 
Y Define Setpoint: what you are trying to achieve 
Y Define Sensing Variable: what it actually is 
Y Define Control Variable what you are actually changing (e.g. power) 
Y Done as a Register to be used elsewhere in TD 
Y Define upper/lower bounds on Control Variable (e.g. max and min 
power) 
Y Generally a good idea to Prevent Integral Windup 
Y Apply output variable to some TD object, usually a heat load 
vY Make sure to output as expression so that changes in process 
variable as SINDA runs are applied correctly 


Heater Edit Form 


Enabled. 


Name: 


Logic Submodel: = MAIN 


Qu gh 
K 
Quow 
Tew Tigh 
Sense Method... Program Qi, i g h 
Use MLI/Insulation nodes if possible 
Apply Heat To: Sense Temperature Fi 
Quad Elem[WFI_OB]:D6FA Top 
Qi ow 
New Tigh 
#lZhe i PAIR OIAR, = fe (Zhe i DI ER ofA eee, 
{Cancel (Help. } 
PID Controller 1 


Enabled... 


Run in Steady State 


Comment: Mid PID Controller for Cryocooler 


BUILD MODEL: Phase Change Material Sig eee * 


¢ Phase change happens over a constant 

temperature 

— PCMs (paraffins) are generally very poor 
conductors of heat 

— PCM assembly is generally a hermetically 
sealed aluminum housing to which a core 
for uniform spreading is bonded 

— Gaps in between core is filled with PCM 

— Core may contact multiple faces for 
improved through conductivity (aka vias) 

— Melting point is material specific. Paraffins 
typically used with Solid to Liquid phase 
change 

ThermalDesktop includes capabilities to 

model PCMs via FUSION function 

Y Define conductivity of material as the core 
(typically 10% or AL1100 for plane of fins, 
and 1% or AL1100 for out of plane) 

Y Define Specific Heat of material as FUSION 

VY Specify Melt Pointing and Heat of Fusion 

Y Majority of mass of Core+PCM is in the 
paraffin for density 


Alkane 
Decane 
n-Undecane 
n-Dodecane 
n-Tridecane 
n-Tetradecane 
n-Pentadecane 
n-Hexadecane 
n-Heptadecane 
n-Octadecane 


Data from http://webbook.nist.gov/chemistry/ 


Comment: 


Formula 
CioH22 
CiiH24 
Ci2H26 
Ci3H28 
CiaH30 
CisH32 
CisH34 
Ci7H36 
CisH3g 


MP (°C) 
-29.6 
-25.6 
-9.6 
-5.4 
BES 
10.0 
18.2 
22,0 
28.2 


Set Col lor... 


Isotropic 


Ise Pressure Scale: 


Aeusion 


(kJ/kg) 
202.3 
142.9 
214.6 
155.6 
224.2 
161.9 
226.9 
165.4 
240.0 


Core material helps to spread heat uniformly 
throughout assembly. Void areas of core are 
filled with PCM. PCM generally stays in 
contact with core 


Liqud-Solid Phase Change 


2202 
Melting Point: 96 
Heat of Fusion: 214629 
Solid-Solid Fusion 
(F Track Solid/Solid2 Fusion 


Solid2 Specific Heat 


Solid-Solid Heat of Fusion: 
Cancel 


Solid Specific Heat 


dékg/C 


c 


i Model Housing as 2D Rectangles, PCM+Core 
i. as Solid with PCM Material assigned. 
~ Contactor between Core/Housing for bond 


BUILD MODEL: Louver Modeling 


¢ Louvers are devices which passively activate as a function of temperature. They are often 
used to reduce the view of a radiator to space when they get colder and increase the view 
as they get warmer to conserve heater power. 

¢ They function when a bi-metallic spring changes its length and spring constant as a 
function of temperature, rotating a set of highly reflective, parallel blades 

¢ These blades change the view factor from an underlying radiator its environment 

¢ Since Radiation couplings are generally for a fixed geometry and properties, the typical 
Monte Carlo Ray Trace does not account for this variability nor does it know the temps... A 

¢ To model this, it is necessary to modify primarily the Radk to Space and environmental 4 
loads in the thermal model. Radiation Model should be run with fully open properties... yw y 

¢ Account for louver in SINDA model (to first order at least)... lip 
Y In OPERATIONS, loop from 1 to NGTOT for each Louver Node. For each Index... | 


Y CALL CONDAT(GETGMOD(Index),GETGNUM(Index), JTEST, KTEST,ZTEST) 
Y Check if NDNAM(ITEST)=SPACE, NDINT(ITEST)=SpaceNode, NDNAM(JTEST)=Louver Sub, NDINT(JTEST)=Louver Node 


Y Check if NDNAM(JTEST)=SPACE, NDINT(JTEST)=SpaceNode, NDNAM(ITEST)=Louver Sub, NDINT(ITEST)=Louver Node 
V If either condition is met, then flag Index as a Louver G and store Cond number for later reference 
Y Store Fully Open conductor value for this G for future modification 

Y In VARIABLES1, for each Louver node (G Index and Fully open conductance should have been determined and stored) 
Y Determine the scale factor for partial closure based on temperature and Louver Performance curves: SCL (1 to >0) 
vY Set G(Louver Index) = G Fully Open * SCL 
v¥ Adjust Heat Load by SCL (If Nadir pointing, Tot Env - Eclipse Env = UV Env, Eclipse Env = IR Env, adjust based on a,s) 


BUILD MODEL: Integrating Delivered Models 


¢ Models delivered from other organization should have already been checked out. Delivered models (PARTs) and assembly 
models (ASSY) often need to be prepared for integration 


RON RRR COR 


~ 


v 


Remove boundary conditions from PART model that represent interfaces: boundary nodes, blocker surfaces. Keep 
conductors/contactors to boundary nodes to be able to connect to non-boundary interface. All that should remain is what needs to be 
imported into next higher assembly. Ensure that Radiation Group names are predefined or do not conflict with those that exist in the 
ASSY file. This should be a clean PART file for import 
Remove previous instance of PART from the ASSY model if necessary or consider if better to keep old PART in ASSY and rename 
submodels with OLD_ prefix to have previous integration template. Save this as a clean ASSY file for integration 
Ensure both PART and ASSY are using the same units 
Use AutoCAD’s INSERT command to insert PART file. Explode and place at correct location and orientation. Could use Copy/Paste if no 
Domain Tag Sets need to be preserved 
If you need to move the entire inserted PART, be careful with any nodes/surfaces that may be attached to assembly. Moving everything 
will move the attached entities twice (relative to assy and then the assy itself) 

Ensure that all PART submodels added have the same number of Nodes, Surfaces, Contactors, Heat Loads, etc. in the 

PART file and the new ASSY file 

Import Optical and Material Properties (Model Browser will show Property Not Found for undefined properties) 

Import Symbols (Model Browser will show Symbol values as /] when referenced but not defined) 

Import Logic Object Manager objects and PART file node correspondence if needed 

Verify that Radiation Groups are correct. Can merge imported PART groups with existing ASSY groups 

Ensure Domain Tag Sets from PART are included in ASSY 

Determine if any symbol or property alias over rides from PART case set are necessary in ASSY case sets 

Reconnect contactors/conductors to ASSY side entities and remove any temporary PART boundary conditions. Merge 

any nodes that may be coincident between PART and ASSY prepared files 

Add Tab to Notes (Utilities-Notes-right Click Tabs) and paste text documentation (text only, no images, tables, etc) 


CHECK MODEL: General Model Checks 


The following list can be used for your own model or models you 


, Submodel -> Submodel A Submodel ... Submodel Z 
received from others: 
. . . 2 Build Geometry X X X 
Y Ensure Material and Optical properties have correct values defined aloo 7 
Y Ensure all surfaces have material and optical properties assigned Nodal Subdivision X 
Merge Nodes X 

Y Verify thicknesses (Model Checks-Color by Property Value-Thickness) Subsea X 
Y Check Mass of Submodels (Model Checks-Calculate Mass) Nodes umbenig a 

Thickness X 
Y Check Free Edges of Elements (Model Checks-Show Free Edges) Material X 
V : ; ; Material Orienters 

Check for Duplicate Nodes (Model Checks-List Duplicate Nodes) Repeae aee 

Y Check for Coincident Nodes (FD/FEM Network-Merge Coincident Nodes) Opt Props 

Insulation 
Y Active Sides for Radiation Group (Model Checks-Display Active Sides) Comment 
Y MLI Assignments (Model Checks-Active Display Preferences) assembly Meche 

IF Contactors 

Y Verify All Contactors make contact (Model Checks-Show Contactor Markers) IF Conductors 

Heat loads 
Vv Verify User Defined nodes are connected (Model Browser) near 
Y Verify heater setpoints and Steady State behavior If building a lot of new geometry, it can be helpful to 

. : : ee have a per submodel checklist to keep track of 

Y Verify proper behavior of assemblies/trackers (Orbit-Display Current Orbit) whether you have fully assigned everything needed 
Y Verify proper associations of Symbols (Mode! Browser) 
Y Verify reasonable values for conductors, heaters, contactors, etc (units check and avoid large couplings) 


CHECK MODEL: Model Browser 34 


r 


¢ The Model Browser is a powerful capability to view the Sst Pecelns-Arabes.Theimslsis (ata | | ©) WARST PreshfLDiscipine AnajsesThermalovg (|) 


ethandineipeaPativommattete bist Edit em Displayes Options Help 

relationships between all the ThermalDesktop objects =~ ors |__| | PRI RAISER Te fen tales (> IL : 
and manipulate what is displayed in AutoCAD Ansys Goup 1| | |¢s om” 
— Many options for listing and controlling display poe L om fa een sieee 
— Shows relationship hierarchy (e.g. node-surface, Conductor-node) eres eae 

¢ Submodel: for model/node organization Aasernbest Teaches peated 

* Non Graphical: Orbits, Properties, Case Sets ea +e 

¢ Analysis Group for Active Side Visualization ee al wexaso no 

* Optical and Thermo Props for property assignments nae eee 

¢ Surfaces/Solids for Conduction Submodels —— Po 

* Contact/Contactor for displaying associated surfaces = i —eomergs 

¢ Assemblies/Trackers for showing assembly hierarchy san os 

* Conductors for user defined conductors warping a . 

¢ Heaters/Heatloads for used defined dissipations a aw 

* Orienters for anisotropic conduction ste 9 epaa 

¢ Fluid objects aa L- a bee an : 

¢ Meshers: TD Mesh definitions act ecto ee 

¢ Mesh Displayers: FEM Mesh Mappers (STOP) a Slee 

* Symbols for showing object dependencies ie a aati 

¢ Groups/Domain Tag Sets for user defined collections sus 


¢ Layers for AutoCAD layers for visibility control ee ke - = 


yy VS SY Sw TN ~ 
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EXECUTE MODEL: Defining Analysis Cases 7 


Define Orbits 
Define Case Set name and File Names 


Define Radiation Tasks (Radiation Groups, Radk, 
HR, Articulating Radks, Orbits) 


Symbol Overrides 

Property Overrides 

Operations Block (TIMEND, Steady, Transient, etc) 
Verify Control Cases 

Define Outputs and Intervals 

Add Comments to document Case Set 


Temporarily Set all Symbols, Aliases, and Property 
DataBases to case set and verify model is as 
intended. Restore back when done verifying... 


Run Case(s) 


Post Process Results 


Heating Rate Case Manager 


Current Heating Rate Case: Op67_AU_Venus 


f 
Orbit: 0p67_AU_Venus 


Basic Orbit | Orientation | Positions | Planetary Data | Solar | Albedo | IR Planetshine | Fast Spin 


Beta Angle: 90} Degrees (Angle between solar vector and plane of orbit) 


{Update | 
Calculated orbital period = 5676.98 sec 


Display Orbit Done | Help | 


|| | Editing 1 Case Set - SIRSE_Decon_to_Sun 


Case Set Manager 
Z) jon | 
All Case Sets efi 
SIRSE_Op et D 
SIRSE_Decon se S 
SIRSE_NonOp Ca 
aes a : 
ee 
l Set Symbols, Aliases, Override Databases to Global ] 


Editing 1 Case Set - SIRSE_Decon_to_Sun 
Radiation Tasks | Calculations | Output | SINDA | Dynamic | Advanced | Props | Symbols | Con 


Optical Properties 


Thermophysical Properties 


| Alias =| 


Override current database 


Radiation Tasks | Calculations | Output | SINDA | Dynamic | Advanced | Props | Symbols | Comments 


Analysis Group Orbit Cale Method Oetens 
SIRSE_Final tk @ Re-use calculated data if valid, 
SIRSE_Final 0p67_AU_Venus br recalculate 


ony n Recalculate data (current 

f i n it | database will be replaced) 

fe S kK D Add rays to database #f possible, 

otherwise recalculate 

{accura cy of current database will 
be refined) 


gadiation 


Always reuse data (no testing 
performed) 


Editing 1 Case Set - SIRSE_Decon_to_Sun 


Radiation Tasks | Calculations | Output | SINDA | Dynamic | Advanced | Props | Symbols | Comments | 


Global List 


Overtide Global 
hrBetaAngle 0 IR_Decon_Power 0 

hrEccen 0 Motor_Dissip 00 095 

hilum 0 VIS_Decon Power 5 0 

hrMean An 360 

| hePeriod 5676.98 

hrPeriodSec 5676.98 

hrPlanetX 4 — 

|hrPlanetY 2.65343e-014 k= 

hrPlanetZ fy 

hrPos 12 

heShadowEntry 4 

hrShadow Ext 4 4 e S 
hrSpeed 7.61261 4 | 
hrSub Solar Mean An 5 e i¢ 
hrSubSolarTrueAn 45 V 
bhrSunX f) 

hrSunY 4 — 

heSunZ 2.65343e-014 

heTime 5676.98 

hrTimeSec 5676.98 

|heTimeShadowEntry 1 

hrTimeShadowEntySec -1 

heTime Shadow Exit 4 

|heTimeShadowExtSec 1 

heTrueAnom 360 

hrVeloctyX a 

hrVeloctyY -2.129922-02 

hrVeloctyZ 1 

Temp_SC_IF 293 


Drive Symbols From Excel 


EXECUTE MODEL: Defining Orbits » 


¢ Orbits are defined to specify the variability of heating based on orbiting a _[Hatemceetense 
celestial object (Planet, Sun, Moon, etc) See at ED 


ee DITL TRAJECTORY Add |__Import | 
¢ Sources typically include: Direct Solar, Planetary Heating, and Albedo GIRO cay BS — 


Cyc6_B90_L2 7-36 S15 Hot —BAS' 


(Solar Reflections off Planet) CO SP BAS 
¢ Orbits can be specified by Altitude, Planetary Data (Size, Mass, eee ze 


Temperature, Fluxes, etc) or by Vector List (Good for small body orbits GeeBSLYE Mca Ba 


Cyc6_B90_L2_Y18_X-15 Cold BAS 
H 


bit manage! z 


Q 


Set Current 


(comet, asteroid)). Orbits internally decompose to vector list for calculations See eS 
¢ Can use Beta = 90 with no Earth/Albedo for Lagrange orbits are 
¢ Orbit Manager allows you to define multiple orbits and specify the current 


LG 
wn 
isivlelelvlelelelelelelelelelelslelelsl=| rn] 


oo as 


Heating Rate Case: CGI_DITL (ee) 


orbit for visualization [Vector ue Planetary Data | Solar_| Albedo | IR Planetshine | Fast Spin 


Ratio of distance from planet center to 
radius 


¢ Orbits are referenced in Radiation Tasks in the Case Set Manager a ee 


0 
111399 =—-0.3211, 0.0000, 0.9470 0.0000, 0.0000,-1.0000 2.000000 


111400 0.2433, 0.0550, 0.9684 0.0000, 0.0000,-1.0000 2.000000 


¢ Orbits can be imported from and exported to other models moe Saeree tae Yee ee ames 


111799 0.1638, 0.1117, 0.9801 0.0000. 0.0000.-1.0000 + 2.000000 
0.1688, 0.982 a 


py oo 111800 0.0833, 0.9821 0.0000, 0.0000, -1.0000 2.000000 
Orbit: Beta_0_LEO Flesting Rate Case orbit 111999 0.0833, 0.1688. 0.9821 0.0000, 0.0000:-1.0000 2.000000 
= a = 3 | 112000 0.0023, 0.2250, 0.9744 0.0000, 0.0000,-1.0000 2.000000 
jasic | Ori Posi Pi IR Pl i F; in | " = r 141999 0.0023, 0.2250. 0.9744 0.0000, 0.0000.-1.0000 2.000000 
Orientation | ositions lanetary Data Solar Albedo IR Planetshine | ast Spin | Keplerian Orbit | Orientation | Positions | Planetary Data | Solar | Albedo | IR Planetshine | Fast Spin | 1 * 03048 09737 0.0000. 0.0000.-1.0000 2.000000 

142199 0.0384, 0.2248, 0.9737 0.0000, 0.0000,-1.0000 2.000000 

Orbit Definition Orbit Shape Definition - Select 2 142200 0.0745, 0.2243, 0. S717 0.0000, 0.0000, -1.0000 2.000000 

— =—_—<— 142399 0.0745, 0.2243, 0.9717 0.0000, 0.0000,-1.0000 2.000000 

Beta Angle: o Degrees (Angle between solar vector and plane of orbit) Orbit Inclination fo] Degrees {[/] Minimum Altitude: 500 km 142400 (0.1103, 0.2236, 0.9684 0.0000, 0.0000,-1.0000 2.000000 
142599 0.1103, 0.2236, 0.9684 0.0000, 0.0000,-1.0000 2.000000 

Attud 00k FA.ct Ascending Mode [80 Degrees EMerimam Aude: (500 hn tg2s89 01481, 0.2228, 09639 0.0000. 0.0000.-10000 2 000000 
sata ™ 192600 0.1461, 0.1491, 0.9780 0.0000. 0.0000:-1.0000 + 2.000000 
Argument of Periapsis: 270 Degrees ¥| Eccentricity: 0 192649 0.1461, 0.1491, 0.9780 0.0000, 0.0000,-1.0000 2.000000 

192650 0.1461, 0.0747, 0.9864 0.0000, 0.0000,-1.0000 2.000000 

Period 5676.98 192699 0.1461, 0.0747, 0.9864 0.0000, 0.0000,-1.0000 2.000000 

Update _ | Date Dependent Right Ascension Definitions — = em : Lyi a — com bom e om Eye 

Calculated orbital period = 5676.98 sec @ User Specified [ Update 192750 0.1461, 0.0747, 0.9864 0.0000, 0.0000,-1.0000 2.000000 


192799 0.1461, -0.0747, 0.9864 0.0000, 0.0000,-1.0000 2.000000 
192800 0.1461, -0.1491, 0.9780 0.0000, 0.0000,-1.0000 2.000000 


Sun 0 Degrees 192849 00.1461, -0.1491, 0.9780 0.0000, 0.0000,-1.0000 2.000000 
. 192850 0.1461, -0.2225, 0.9639 0.0000, 0.0000.-1.0000 2.000000 
it P, 242900 0.1461, 0.2225, 2.9639 0.0000, 0.0000,-1.0000 2.000000 
( rime Meridian 0 Degrees 


C O pit : 
pas! a —— xeplerian or or List 


Calculated Beta Angle: 0.000000e+000 


———— (Gee) Gate) | | Ss 


LN ANNAN 


EXECUTE MODEL: 


Open Orbit Manager and select Add... 
Define Type of orbit and Specify Name 
Enter Basic orbit information (e.g. Alt., Beta Angle, Incl., RAAN, Vector List, etc) 
Specify orientation of SC wrt celestial objects (e.g. Sun/Nadir/Zenith Pointing) 
Specify Euler Rotations to orient SC. Can use hrMeanAnom to do slow spin... 
Define number of orbit positions (default is 12: 15 positions with beginning and 


end + Eclipse Entry and Exit) 


~ 


Specify environmental heating (Solar, Albedo, IR Planetshine) Albedo and 


Planetshine can be Lat/Long dependent 


Orbit: Beta_O0_ LEO 


Basic Orbt | Orientation | Positions | Planetary Data | Solar | Albedo | IR Planetshine | Fast Spin 


7 


Orbit: Beta_O_LEO 
Beta Angle: 0] Degrees 


Atitude: 500 km Pointing 


|_Update _] 
Calculated orbital period 


Star 


Right Ascension: 0 Degrees 


Declination: 0 Degrees 


Velocity vector 


Orientation Override 


Align to Celestial Coordinate System 


Basic Orbit} Orientation | Positions | Planetary Data | Solar_[ Albedo [ IR Pianetshine | Fast Spin | 


Additional Constraint 


ne: (Sa 


Nadir 


Star 


Right Ascension: 


Declination 
®@ Velocity vector 


Additional Rotations 


aa| (90 
[x | -360 
[Zz] 0 


0 q 
0 q 
Deg 
Deg 
Degr 


ta 


Defining Orbits 
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a 
Create New External Heating Environment 


New Heating Case Name: 


Orbit: Beta_0_LEO 


Basic Orbit | Orientation | Positions | Planetary Data | Solar_| Albedo | IR Planetshine | Fast Spin 


Type 
® Basic Orbit 
Keplerian Orbit 
Heating Rate Case Mar : e é A, 
ane Planetary Latitude/Longitude/Altitude List 
Current Heating Rate Case: Beta_O_LEO : : A 
Vrms Orbital Sun/Planet/Radius Vector List 

CGI_DITL TRAJECTO ; : 2 
CGI_DITL_Visual BASIC Free Molecular Heating Velocity Vector List 
Cycb_1Day_Orbit BASIC 


Use Positions: 


Shadow Crossings (Degrees) 
Start: 0 Degrees 
¥| Automatically Include: 
Entry: 111.984582 
End 360 Degrees a 


Ext: 248.015418 
Position Angles 
@) True Anomaly 


fn Mean Anomaly 
[ Update List. | 


Edit list items directly, use <Enter> to insert entry 


Edit Planetshine vs. Time Table... 


Use Planetshine vs. Latitude/Longitude 


Edit Planetshine vs. Lat/Long Table... 


®@ Planet Coordinate System 


Subsolar Coordinate System 


Cyc6_B90_L2_Y-18.X15_Cold BASIC 
Cyc6_B90_L2_Y-36_X-15 Cold BASIC Seay - 
Cre BOTLEY IEKIEHot BASIC Free Molecular Heating with Reference Orbit 
Cyc6_B90_L2_Y-36_X0_Cold BASIC 
Cyc6_B90_L2_Y-36_X0_Hot BASIC H i i i i 
CyeR BONY IES Cold BASIC Celestial Coordinate System Location/Orientation 
CycB_B9O_L2_Y-36_X15_Hot BASIC. 
Cyc6_B90_L2_Y0_X-15_Cold BASIC 
CycB_B9O_L2_YO_X-15_Hot BASIC 
Cyc6_B90_L2_YO_X0_Cold BASIC ) 
Cyc6_B90_L2_YO_X0_Hat BASIC | ...0K...] 
Cyo6_B90_L2YO_X15 Cold BASIC uals | Cancels) hele 
Cyc6_B90_L2_Y0_X15 Hot BASIC 
Cyc6_B90_L2_Y18_Xx-15 Cold BASIC 
Cyc6_B90_L2-Y18_X-15 Hot BASIC 
Cyc6_B90_L2_Y18_X0_Cold BASIC 
Cyc6_B90_L2_Y18_X0_Hot BASIC 
CyeB_BSOLZ VIE XTS Cold BASIC : 
Display Orbit Done | l Help 
==] 
Orbit: Beta_0_LEO 
Basic Orbit | Orientation | Positions | Planetary Data Solar Albedo | !R Planetshine | Fast Spin 
Options 
Input Mode: 
Value: 0.000243 W/mm*2 Temperature (Black Body) 
0.0002215 @ Flux 
Use Planetshine vs Time Planetshine Coordinate System 


EXECUTE MODEL: Defining Radiation Tasks ‘a 


Radiation Tasks are a subset of a Case Set. The results of these 
Tasks serve as inputs to the thermal model and generally include 
internal radiation within an enclosure, external radiation, and 
environmental heating. Radiation can be for moving geometry. 
Each radiation task generates a file that is included along with the 
Cond/Cap calcs to form the SINDA model for temperature solution 
Results from Radiation Tasks may be used by multiple case sets. 
TD is generally fairly smart at evaluating if anything has changed 
since the last time a radiation task was run and determining if the 
data is still valid. User can over ride this... 

Each radiation task has its own control parameters including 

Max Rays to fire, Acceptable Error Criteria, Output File and 
Submodel, Oct Cell Subdivision, etc 

Fast Spin Capability also exists to vary a spin angle multiple 

times within one orbit position. Averaging these results 
(Radk/HR) together over multiple spin positions in one orbit 
position can simulate a fast spinning object (Scan Mirror, 
Reflector, etc). To effectively utilize this, an assembly should be 
created with the spinning geometry attached and the rotation 
based on the symbol representing the spin angle. Use 0 Error! 


Editing 1 Case Set - WFIRST_Launch_Orbit_BO 
Radiation Tasks | Calculations | Output | SINDA | Dynamic | Advanced | Props | Symbols | Comments 

Analysis Group Orbit Cale Method Options 

CORONA rk © Re-use maar data if valid, 

INT_IC ik oe 

INT_PM rk 

INT_SM rk = ean e data (current 

INT_WFI_IFM ik jatabase will be replaced) 

INT_AMS rk 

WFIRST_Extemal rk Add ic dhsishican paca: 

WFIRST_Extemal Beta_0_LEO hr erases serene agi 
{accuracy of current ~——— will 
be refined) 
Always reuse data (no testing 

Radiation Analysis Data [<3 performed) 


Job | Control. | Advanced Control | Radk Output | Spin _ | Overap 


Calculation Type 
@ Radks: 


Radiation Analysis Data 
Heating Rates 
Job Control | Advanced Control | Radk Output | Spin | Overlap 
Atticulating Radks | ; 
Default 
Free Molecular Conduction | 
Atticulating Free Molecular Conduction | Rays Per Node: 35000 
View Factors Weighted Error: y ee 
. . ° Delete | 
Atticulating View Factors | Rays Before Initial Error Check: 1000 
Energy Cutoff Fraction: 0.01 
Analysis Group WFIRST_Extemal 
Orbit: Beta_0_LEO | Heating Rate Sources 
Calculation Method 
@ Monte Carlo 
| 
Progressive Radiosity Nodes Positions 
@ All @ All 
List List 
Add to Database Name: None 
Radiation Analysis Data [nese] 


Job | Control | Advanced Control | Radk Output | Spin | Overlap 


¥| Do Fast Spin Calculations 


Symbol Name Roll] 


@ Use Equal Incrments 


Starting Value 0 
Stop Value: 360 
Number of Increments: 12 
Inclusive Start and Stop Values 
User Defined List 


Y Open Case Set, Go to Radiation Tasks Tab, Click Add... 


Y Define Type of Task (radk, Heating Rates), Analysis Group, and Orbit if needed 
Y Add to Database Name if you don’t want to overwrite similar 
— Name of Database is RadGroup-Orbit-PropFile.rch or RadGroup-PropFile.rck 


Y Specify Max Rays, Error, and Energy Cutoff under Control! 


Y Specify Levels and Objects on Advanced Control (Objects > 8) 
Y Set output filename and submodel...DO NOT USE DEFAULTS!! 
Y Set Bij Cutoff and Sum and disable Radks to space if internal radiation 


EXECUTE MODEL: Defining Radiation Tasks 
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Editing 1 Case Set - WFIRST_Launch_Orbit_B0 


Analysis Group 
CORONA 

INTIC 

INT_PM 

INT_SM 
INT_WFI_IFM 
INT_AMS 
WFIRST_Extemal 
WFIRST_Extemal 


Orbit 


Beta_0_LEO 


Cale Method 


THereernee 


Radiation Tasks | Calculations | Output | SINDA | Dynamic | Advanced | Props | Symbols | Comments 


Options 


@ Reuse calculated data f valid, 
jerwise recalculate 


Recalculate data (current 
database will be replaced) 


Add rays to database if possible, 
otherwise recalculate 

(accuracy of current database will 
be refined) 


Always reuse data (no testing 


= 


= 
Radiation Analysis Data 


Calculation Type 
© fa 


Free Molecular Conduction 


© View Factors 


Articulating View Factors 


Orbit: 


Calculation Method 


@ Monte Carlo 


Progressive Radiosity 


Fl 


Radiation Analysis Data 


Job | Contro| Advanced Cortrl| Rack Output | Spin _ | Overiap 
i 


performed) 


rr 
Radiation Analysis Data 


Job | Control | Advanced Contral | Hestrate Output [Spin _| 


Job | Control | Advanced Control | Ras 


Articulating Free Molecular Condu: 


Analysis Group: BASE 


Apply Reciprocity To View F 


Add to Database Name: [None 


Radiation Analysis Data 


Job _| Control | Advanced Control | Radk Output | Spin _ | Overap | 


Default 
Rays Per Node: 
Weighted Error: 
Rays Before Initial Error Check: 
Energy Cutoff Fraction: 


Heating Rate Sources 


35000 


Positions 
@ All 


List 


Oct Cells 


[¥] Use Octtree to accelerate calculations 
Max oct-ree subdivisions: 12 


Max surfaces per cell 12 


Random Number Seed Control 
@) Use unique random number seed at start of calculations 
Use same random number seed sequence at start of calculations 


Use same random number seed sequence at start of every node 


Nodalization Schemes 


@) General Specific Top/Bottom Specific 


Radk Calculation Spectrum - Used for Modeling Lamps 
© Infrared Solar 


Wavelength Dependent Properties 


¥] Enable Wavelength Dependent Calcs Edit... 


Radiation Analysis Data 


[Jeb | Control | Advanced Control | Rack Output [Spin _| Overlap 


| Generate SINDA/ECUINT input after calculations: 


Radk Output Filename: 

Form Factor Output Filename 
Output Submodel 

Initial Conductor ID: 

Space Node Submodel ID: 


Space Node Temperature: 


List summary if “kept is off by more than 


Which Radks Output To SINDA 


Bij/Fij Cutoff Factor: 


¥| Output Bijs Until BijiSum >= 


| Output 1-BijSum to ERN Node: 


[©] Disable Radks To Space 


| Output as Heating Rates: 


WFIRST_PreSRR_Extemal.k 
WFIRST_PreSRR_Extemal.dat 


WFIRST_RDK 


Output Filename: WFIRST_PreSRR_Y18_X-15_Hot.hr 
Output Submodel: WFIRST_HR 

S/F Starting Aray ID: 1 

Output Format: LOADQ 


[¥] Combine SAP arrays into a single array 


¥] Output tracker data to file; WFIRST_PreSRR_YO_X0_Hot.dat 
"| Output HR Symbols to SINDA 


[ Advanced Logic... } 


Multiple files will be produced 
"A" will be appended for ARRAY DATA 
"L" will be appended for VARIABLES 0 Logic 


] Output as fluxes 
Sources 
¥| Solar Planetshine | Albedo 
Diffuse Sky Solar | Diffuse Sky IR 


] Diffuse Sky Albedo 


EXECUTE MODEL: Radiation Calculations 


¢ Monte Carlo Ray Trace ° Oct Cells 

— Begin with Surface 1 and find a random point and 6=0.5 — Find bounding box around 
direction to “fire” the ray with 100% energy —jyo% aN. 100% entire model 

— Determine if the ray intersects any surface so, “""\-*sy_— = Divide by midplanes in all three 

— Oct Cells subdivide model to reduce number of 302»~ oe 0.7 directions to form 8 smaller 
intersection tests boxes 

— If ray intersects a surface, determine how much energy is — Determine #surfaces in each of , 
absorbed by the surface (a,6) the 8 cells iy 

— If some energy remains, determine if the energy is — If #surfaces is > specified, then subdivide oN 
reflected specularly (Angle of Incidence = Angle of cell again. Continue until max #subdivisions / 
Reflection) or diffusely (Random direction) or max #surfaces is reached 

— Continue propagating ray until energy is below extinction — Consider 8 surfaces meeting at point (i.e. 
threshold at which point it is either completely absorbed ribs). Lower #surfaces can never be met... 
or completely reflected ° Error Calculations 

— Once extinguished, select new random point and — Statistically, the error can be estimated with 90% 
direction for next ray confidence by: 

— Once finished with Surface 1, move to next surface Werk. = 1.65 1-8 i yx 100 

¢ How Many Rays? i NN ays Bi 

— Typically start with 35000. (2.65/1.76% error for 0.1/0.2 B,) | — Weighting the Error for every node yields a single term 

— Is the run time tolerable? per surface to determine if sufficient rays have been 

— Spending time calculating small couplings that are eliminated? __ fired for that surface 

. n 
— Does it affect temperatures? By how much? Weined tae dj=0(B; * Error) 


In the end it comes down to run time vs. accuracy... jeo(Bi,) 


VERIFY RESULTS 


Verify Results refers specifically to the process of confirming that the model and predictions are consistent with 
the physics of the design. 

Model correlation is a measure of how accurately the model predicts actual measured values under the same 
conditions as tested...this happens well after a model is built, much later in a program 

The first thing to check are temperatures. A contour plot is useful for this as it will highlight any gradients 

If the contours look reasonable (what should be cold is cold, what should be hot is hot), the next thing to 
compare is the temperatures against design limits. XY plots are good to see the local max and min values over 
an orbit, but it is hard to do this for every node of importance. Comparison of the hottest temperature and 
coldest temperature of a component (or group of node) over the orbit to the design limit is best accomplished 
using a spreadsheet and is often represented in a tabular form. 

XY plots are also useful to display the nodal stability and compare against any stability requirements 

For heaters, it is important to confirm that the heater power is sufficient. No more than 70% of the available 
heater power should be used per GSFC GOLD rules. The required amount of power to maintain a component 
within limits does not depend on heater size, but duty cycle does... 

A variety of external tools exist for post processing. TD offers XY, Contour Plots, Animated Contour Plots and 
methods for: querying Max/Mins, writing output text files, and calculating sinks and heat flows post-solution 


Lastly, to really understand a design, there is no substitute for investigating heat flows. Heat flows identify 
areas of a design where better coupling is needed, better isolation is needed, or unexpected heat inputs or 
leaks are present. This takes a fair bit of time to set up, but the understanding of a model/design is far better 


after having gone through the effort. 


